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Abstract. 

The search for signatures of transient, unmodelled gravitational-wave (GW) 
bursts in the data of ground-based interfcrometric detectors typically uses 'excess- 
power' search methods. One of the most challenging problems in the burst-data- 
analysis is to distinguish between actual GW bursts and spurious noise transients 
that trigger the detection algorithms. In this paper, we present a unique and 
robust strategy to 'veto' the instrumental glitches. This method makes use of the 
phcnomcnological understanding of the coupling of different detector sub-systems 
to the main detector output. The main idea behind this method is that the 
noise at the detector output (channel H) can be projected into two orthogonal 
directions in the Fourier space - along, and orthogonal to, the direction in which 
the noise in an instrumental channel X would couple into H . If a noise transient in 
the detector output originates from channel X, it leaves the statistics of the noise- 
component of H orthogonal to X unchanged, which can be verified by a statistical 
hypothesis testing. This strategy is demonstrated by doing software injections in 
simulated Gaussian noise. We also formulate a less-rigorous, but computationally 
inexpensive alternative to the above method. Here, the parameters of the triggers 
in channel X are compared to the parameters of the triggers in channel H to see 
whether a trigger in channel H can be 'explained' by a trigger in channel X and 
the measured transfer function. 



PACS numbers: 95.55.Ym, 04.80.Nn 

1. Introduction 

A world-wide network of gravitational-wave (GW) detectors consisting of ground- 
based interferometers and resonant-bars jSlElCl has now started looking for 
signatures of GWs expected to be coming from astrophysical and cosmological sources. 
Among the most promising astrophysical sources of GWs for these ground-based 
detectors are the transient, unmodelled astrophysical phenomena like core-collapse 
supernovae, Gamma-ray bursts and black hole/neutron star mergers. Although, in 
recent years, numerical relativity has made a tremendous progress in predicting the 
gravitational waveforms from compact binary mergers and gravitational collapses (see, 
for e.g., [8] for some recent work), these predictions are not yet robust enough so as 
to allow the data-analysts to use techniques like matched filtering, which relics on 
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the accurate models of waveforms. Most of the algorithms currently used in burst- 
searches are time- frequency detection algorithms which look for short-lived excitations 
of power in time- frequency maps constructed from the time-series data ^lEHE]- Such 
detection algorithms implemented in the data analysis pipelines are usually called 
event trigger generators (ETGs). 

Since current interferometric detectors are highly complex instruments, the 
detector output will contain a large number of spurious noise transients which trigger 
the ETGs. One of the main challenges in the burst-data-analysis is to distinguish these 
instrumental 'glitches' from actual GW bursts. One way of approaching the problem 
is to use more than one detector (mostly located at different parts of the world) and 
look for coincident burst signals in these data streams. While coincident instrumental 
bursts in multiple detectors are highly improbable, long data-taking runs (typically 
several months long) with multiple detectors can produce a large number of random 
coincidences. All of the detector groups are thus developing various techniques to 
'veto' the spurious instrumental bursts at the detector site itself in order to reduce 
the list of candidate coincident triggers (see |12l 1131 IT^| for some recent work in this 
direction) . 

In this paper we present a robust veto strategy which makes use of a 
phcnomenological understanding of the coupling of various detector sub-systems with 
the main detector output. This provides us the possibility of vetoing the triggers with 
a very high confidence. Although this method was developed for the GEO 600 [3] 
detector, this can be applied to the data from any interferometric GW detector. 

The main idea behind this method is that the noise measured at the detector 
output can be decomposed into two orthogonal components. If the physical coupling 
between a sub-system X and the detector output is measured, this information can 
be used to transfer, or map, the noise in X to the detector output, and hence the 
component of H orthogonal to the 'mapped' X can also be calculated. If a noise 
transient originating from the sub-system X appears in the detector output, it leaves 
the statistics of the noise in the direction orthogonal to the mapped X unchanged. In 
SectionsElandElwe formulate and demonstrate a strategy to veto such noise transients 
using the transfer function from X to H. In Section^ we develop an alternative veto 
strategy using the parameters of the burst triggers estimated by the ETG. 



2. Vetoes using known instrumental couplings 

Let x{t) and h(t) denote the time-series data measured at the input and output of a 
linear, time-invariant system. The input and output of the system are related by the 
transfer function Txh(s) of the system, defined as 15 

Txn(S) = "TTTTj (1) 

X(s) 

where H(s) and X(s) represent the Laplace transforms of x(t) and h(t), respectively, 
and the complex variable s represents a point in the Laplace space. Although the 
transfer function Txh(s) is formally defined in the Laplace space, for the purpose of 
this paper, it is easier to work in the Fourier domain. If X(f) and H(f) are the 
Fourier transforms of x(t) and h(t), respectively, the equivalent relation in the Fourier 
domain is given by 

Tx„(/) = §^, (2) 
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Txh(/) = ^, (3) 

where Pxh is the cross-power spectral density of x(t) and h(t), and Pxx is the power 
spectral density of x(t). In the formal language of signal processing, Txh(/) is 
described by the frequency response (magnitude of Txh(/)) an d phase-shift (phase 
of Ixh(/))- But we will conveniently think of Txh(/) and Txh(s) as two different 
representations of the transfer function. 

We will refer to the measurement points for time-series data within the detector 
as 'channels' and assume that they are continuously recorded. Let x(t) and h(t) 
denote the time-series data in channel X (which, presumably records the noise from 
a detector sub-system) and channel H (the main detector output), respectively. The 
transfer function from X to H can be measured by injecting some noise in X and 
measuring X(f) and H(f) simultaneously ^Hl- This is done in such a way that the 
injected noise from X completely dominates channel and the contributions from 
other noise sources are negligible. The measured transfer function represents our 
phcnomenological understanding of the physical coupling of a detector sub-system to 
the main detector output. If we assume that the coupling of noise between channel 
X and H is linear and the transfer function is time-invariant, the Fourier transform 
of noise measured in channel X at any time can be transferred to channel H, by 

X'(f)= T XH (/)X(/). (4) 

X(f) and H(f) can be thought of as components of the vectors X and H, 
defined in two different infinite-dimensional Hilbcrt spaces spanned by the Fourier 
basis functions. In that sense, Eq.Q is equivalent to mapping X into the space of H. 
Then, the component of H that is orthogonal to X' can be found by a Gram-Schmidt 
orthogonalization ^Zj 

£ = H-proj x ,H, (5) 

where we define the projection operator by 

proj a v = -pr-^r u . (6) 
(u, u) 

In the above expression (v,u) denotes the inner product between the vectors v and 
u: 

fCX. 

(v,u) = 2/ «(/)u(/)*d/. (7) 
Jo 

All of the continuous quantities like H(f) and X(f) can be recast in terms of their 
discrete counterparts like Hk and Xk , where the index k represents the fcth frequency 
bin of the discrete Fourier transform (DFT). Now H and X can be thought of as 
vectors defined in two iV-dimensional Hilbert spaces, where N is the number of data 
samples used to compute the DFTs. 

The noise in different channels is generated by different random processes and 
has specific statistical and spectral properties. If a non-stationarity (e.g., a glitch) 
originates in channel X, it will change the statistical properties of that segment of data 
in channel X, and hence, in channel H . But the statistical properties of 8 constructed 
from this segment will remain unaffected. On the other hand, if the non-stationarity 
does not originate in X, it will change the statistical properties of 8 constructed from 
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this segment. This can be verified by a statistical hypothesis testing. If S constructed 
from the segment of data containing the non-stationarity is statistically similar to 
that constructed from the neighboring segments, this means that the non-stationarity 
originates in an instrumental sub-system (channel X) and we veto the trigger. If this 
is not true, this means that the non-stationarity does not originate in channel X, and 
we keep the trigger. This exercise can be repeated with all the known noise sources 
(all the measured channels that are known to couple to the detector output). 

2.1. Test statistic 

If the components of H are generated by zero-mean Gaussian processes, the real 
and imaginary parts of 5k in each frequency bin will be distributed according to a 
Gaussian distribution of mean zero and variance er|. Following we compute the 
'excess-power' statistic from 5: 



It can be shown that e will follow a x 2 distribution of 2M degrees of freedom in 
the case of a non-windowed DFT. But in the case of a windowed DFT, are not 
independent x 2 variables, and hence e will not follow a x 2 distribution |18| . But to 
a very good approximation, e will follow a Gamma distribution with scale parameter 
a and shape parameter (3. These parameters are related to the mean and variance of 
the distribution of e by 



In order to estimate the parameters of the Gamma distribution, we generate a 
population of e from stationary data surrounding the burst (using the same DFT- 
length and bandwidth). From that population, /i e and of can be estimated, and 
hence a and fi. 

If the computed e (from the segment of data containing the burst) is less than a 
threshold, we veto the trigger. The threshold r giving a rejection probability of ijj can 
be found by 



Jo 

where r(x; a, (3) is the probability density of the Gamma distribution with parameters 
a and (3. 

3. Implementation 

Two sets of burst triggers are generated by running an ETG on channels X and H. 
We take a set of triggers that are coincident in channels X and H, allowing a liberal 
time-window for coincidence. The data is divided into N number of segments, each 
of length L. The test statistic, e, is computed from the segment of data containing 
the burst. It is well-known that the maximum signal-to-noise ratio (SNR) for the 
'excess power' statistic is achieved when the time-frequency volume used to compute 
the statistic is equal to the actual time-frequency volume of the signal [§]. Since the 
duration and bandwidth of the burst is estimated by the ETG itself, this information 




(8) 




(9) 




(10) 



Instrumental Vetoes 



5 




is used to decide on the length (L) of the data-segment used to compute 5k and 
the bandwidth over which the integration is carried out in Eq.(JHJ. Consequently, 
the frequency resolution of the DFT used in the analysis is in general different for 
each trigger, and hence, so arc the dimension of the vectors H and X. It is then 
required that the discrete transfer function vector should also have the same dimension. 
So we store a high-resolution transfer function and interpolate it to the required 
frequency resolution. It was found that the analysis can be sensitive to the errors 
in the interpolation, since the interpolation can smear out the detailed features in the 
transfer function. Since, the lower the frequency resolution the higher are the errors, 
we set up a minimum frequency resolution of 16 Hz for the analysis. In order to 
achieve this, the minimum length of the data segment used to compute the DFT is 
set to be 1/16 s ~ 60 ms. 

The parameters of the Gamma distribution are estimated from segments of data 
neighboring the one containing the burst, but excluding that segment. The trigger is 
vetoed if e < r, where the threshold r giving a particular rejection probability ip is 
calculated using Eq. l|10[l . 

3.1. Software injections 

Let us define some terminology. The rejection probability is the probability of a trigger 
originating from channel X being vetoed using the method described above. The 
veto efficiency is the fraction of such triggers that are actually vetoed. In order to 
demonstrate that the obtained veto efficiency is in good agreement with the rejection 
probability, we simulate a population of bursts in channel X and H such that these 
are consistent with the transfer function from channel X to H . We then try to veto 
these triggers after choosing a specific rejection probability. If all of our assumptions 
are valid, the fraction of vetoed triggers among this population should be equal to the 
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Figure 2. Amplitude spectral densities of X, X', H and 8 in the presence of a 
sine-Gaussian injection. It can be seen that 8 contains no trace of the injected 
signal. The injected sine-Gaussian has a central frequency of 1715 Hz and SNR 
of 140 in channel X. 

rejection probability. 

We generate a data stream of Gaussian white-noise and 'inject' Gaussian- 
modulated sinusoidal waveforms in to it; this forms our channel X . This data stream is 
filtered using a time-domain filter and some extra noise is added to it. This constitutes 
our channel H (the 'extra' noise being the component of H orthogonal to X'). The 
frequency response of the filter is the transfer function from X to H. The transfer 
function used in this simulation is shown in Figure ^ which is quite similar to one 
particular transfer function measured in GEO 600. The injected sine-Gaussians are 
of the form: 

s(t) = s ISS (^pj ' sin [27r/o(t - to)] cxp [-(* - hf/r 2 } , (11) 

where /o is the central frequency of the waveforms and to is the time corresponding 
to the peak amplitude. We setup the envelope width as r = 2//o. The corresponding 
quality factor is Q = v2wfoT = 8.9 and bandwidth is A/ = fo/Q ~ 0.1/n. The 
quantity s rss is the root-sum-squared (RSS) amplitude: 

1/2 

— s rss- (1^) 

The s rss is varied so that the SNR \ of the injections ranges from ~ 6 to ~ 500 in 
channel X, and the central frequency is randomly chosen from the interval (432 Hz, 
3008Hz). 

As an illustration, the amplitude spectral densities of X (data vector in channel 
X), X' (data vector in X, mapped to H), H (data vector in channel H) and 8 

\ The SNR, p, of a burst trigger is defined by p 2 = 2 / ( \s(f) | 2 /|n(/) | 2 ) d/ where s(f) and h(f) 
represent the Fourier transforms of the signal and the underlying noise, respectively, and the 
integration is carried out over the bandwidth (positive frequencies) of the burst waveform. 
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Figure 3. [Left]: Estimated veto efficiency plotted against the rejection 
probability. [Right]: Estimated false- veto probability plotted against the rejection 
probability. 

(component of H that is orthogonal to X') in the presence of a particular sine-Gaussian 
injection are shown in Figure The injected sine-Gaussian has a central frequency 
of 1715 Hz and SNR of 140 in channel X. As expected, S contains no trace of the 
injected signal. 

The fraction of vetoed events among the injections is plotted against the rejection 
probability corresponding to the chosen threshold in Figure (left). It can be seen 
that the obtained veto efficiency is in very good agreement with the chosen rejection 
probability. 

There is a non-zero probability for an instrumental burst which does not originate 
in A, or a true GW burst, to be vetoed using this method. But, since the probability 
densities of neither GW bursts nor instrumental bursts are known a priori, there is 
no rigorous method for estimating this 'false-rejcction'/'falsc-vcto' probability. As 
a plausible estimation, we inject two populations of sine-Gaussian waveforms with 
random parameters into two data streams of white-noise (so that the waveforms in 
channel X and H are inconsistent with the transfer function from X to H). We then 
try to veto these triggers using the transfer function from X to H . The estimated 
'false-veto fraction' is plotted against the chosen rejection probability in Figure 01 
(right). This suggests that veto efficiencies of > 92 % can be achieved with a false- 
veto probability of < 1%. 

4. An alternative method: 'trigger mapping' 

Although the above described method is rigorous and makes use of the complete 
information contained in the data, the method can be computationally expensive 
because it involves reprocessing of the time series data from the two channels. It 
may be noted that the whole data stream is processed by the ETG in the first place 
and a condensed form of the information about the burst waveform is stored (which 
is often known as 'metadata'). In this section, we develop a strategy to veto the 
spurious triggers in channel H by comparing them with the set of triggers in channel 
A, making use of the transfer function from X to H as well as the information 
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Figure 4. A schematic diagram of the parameter space of burst triggers. The 
cross represents a trigger in channel H and the circle represents a trigger in 
channel X, mapped to channel H. A distance vector w is drawn from one point 
to another. The consistency volume chosen to veto the trigger is also shown. 



extracted by the ETC Although this method is not as rigorous as the previous one, 
the advantage is that this does not require the reprocessing of time-series data and 
hence is computationally inexpensive. 

Let {Ex} and {Eh} denote the set of burst triggers in channel X and channel H, 
respectively. Let us assume that each event, E, is parametrized by its central frequency 
, amplitude E a and time-of-occurrence E l . It is useful to think of E as a point in a 
three-dimensional parameter space with coordinates (E a , , E 1 ). Using the transfer 
function from X to H, we can predict how a certain event Ex would appear in H. 
In other words, we map the event Ex to the space of Eh, making use of the transfer 
function from X to H. In order to veto an event Eh in channel H, we check whether 
any of the 'mapped' Ex triggers are consistent with Eh in time-of-occurrence, central 
frequency and amplitude. Indeed, the precise definitions of these parameters depend 
upon the ETG, and hence we make use of these definitions in order to map the burst 
triggers from one channel to the other. 

4-.1. Mapping the burst triggers 

The burst detection algorithm called HACR |19l 120] is used to generate the burst 
triggers. HACR makes a time- frequency map of the data and identifies time-frequency 
pixels containing excess power which are statistically unlikely to be associated with the 
underlying noise. HACR then proceeds to cluster the neighboring pixels containing 
excess power to form an 'event'. The central frequency and time-of-occurrence of the 
burst-event are estimated by a weighted averaging of the pixel coordinates. This is 
equivalent to the calculation of the center-of-mass of an extended object where the 
signal power in a pixel serves as the 'mass' term. HACR also estimates the total power 
contained in all the pixels belonging to a particular event, and the peak-power of the 
event. The square root of the total power is taken as the characteristic amplitude E^ 
of the event. 

The ETG docs not reproduce the complete physical properties of a burst 
waveform. Instead, the ETG tries to parametrize the underlying waveform using 
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Figure 5. Distribution of HACR errors (after subtracting the mean and 
normalizing by the standard deviation) in estimating the parameters of the 
injected sine-Gaussian waveforms. Also plotted is the probability density of 
Normal distribution with mean and variance 1. 



a set of quantities like the central frequency, amplitude, bandwidth, duration, etc. 
Considering the fact that we are mostly looking at short-lived, band-limited bursts, 
we approximate the power spectrum of the underlying burst waveform in channel X 
to a Gaussian function. For example, the power spectrum of the waveform associated 
with a trigger Ex is approximated to a Gaussian function G(f) such that 

/ /2 G(/)d/=(^) 2 , (13) 
Jfi 

where the limits of integration arc defined by bandwidth _E X of the burst, i.e. 



h = E{- E^/2 and / 2 



Ei + E b x /2. 



Since the peak-power of the burst is also 



estimated by the ETG, the 'spread' of the Gaussian function can be calculated by 
solving Ea. l|13|) . The power spectrum is 'deformed' by the transfer function Txn(f) 
when the glitch makes its way to channel H] this we denote by 



G(/)=G(/)|T XH (/)| 5 



(14) 



Given the transfer function Txh(/), the burst triggers in channel Axan be mapped to 
channel H in the following way: 

1/2 



Eg 



G(/)d/ 



iW — -C/y 



f f f *G(f)\(f)df 
^ 2 G(/)d/ 



E- 



fi 



f£G{f)fdf 
// 2 G(/)d/ 



(15) 
(16) 
(17) 



In the above expression, A(/) is the frequency-normalized phase-delay (time-lag) of 
the transfer function Txa(f)- i-C, 



A(/) = ^0(Txh(/)). 
where 0(.) denotes the phase of a complex quantity. 



(18) 
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Figure 6. Standard deviation of the distribution of the HACR errors (in 
estimating the parameters of the injected sine-Gaussian waveforms) plotted as 
a function of the SNR of the triggers. Also shown are the polynomial fits to the 
data. 



4-. 2. Identifying consistent events 

In order to veto a trigger, Eh, in channel H, we check whether any of the 'mapped' 
triggers (E x ) from channel X, are 'sufficiently close' to it. To be explicit, we define a 
vector, w, connecting the two points in the parameter space of H triggers, 

w = E H -E x , (19) 

and require that it has a sufficiently small 'length'. This length is assigned to the 
vector by calculating the fractional volume enclosed by a three-dimensional Gaussian 
envelope of width cr(w). This is explained below. 

Let AE denote the errors in the ETG in the estimation of parameters associated 
with the event E. In the absence of any systematic biases, the errors AE can 
be assumed to be drawn from multivariate Normal distributions of zero mean 
and standard deviation cr(AE), where the standard deviation is an exponentially 
decreasing function of the SNR. er(AE) can be estimated by injecting known waveforms 
in to the data and comparing the trigger-parameters estimated by the ETG to the 
actual parameters of the injected waveforms. 

Once the errors AEx associated with Ex are estimated, they can be mapped 
to the space of H triggers using Eq. i|15ll7|) by a linear approximation of the error 
propagation |21j . Then the components of w will be distributed according to Normal 
distributions of zero mean and the following variance: 

a 2 (w) =( T 2 (AE H )+a 2 (AE x ), (20) 

where AE X denotes the errors AEx, mapped to the space of H triggers. Thus the 
joint probability density of the vector w comprising of the three random variables 
(w a ,w* ,u>*) is given by the three-dimensional Gaussian function 

■^( w ) = (o ^3/2 i fw a\ i t\ CX P (~^ wTc ~ lw ) > ( 21 ) 
(2ir)- i '' : a(wJ ) a{w a ) a('ur) \ 2 J 

where w T denotes the transpose of w. Assuming that the errors are uncorrelated, we 
write the covariance matrix C as, 

a 2 (w f ) 
a 2 {w a ) 
ct 2 (u>*) 



(22) 
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Figure 7. The RSS amplitude, E a , of the burst-events plotted against the central 
frequency, E?. Shown in the Figure are three population of events: events in 
channel X (Ex); channel X events mapped to channel H (E x ); and events in 
channel H (Eh). 



This enables us to set a threshold for w for vetoing an event Eh- In order to veto an 
event, we require that 

|w| < t. (23) 

The components (r a ,r*,r^) of the 'threshold vector' r are related to the rejection 
probability ip by 

fat 

ip = 8 / / / f(w)dw f dw a dw t . (24) 
Jo Jo Jo 

It can be seen that, by choosing a particular threshold vector x, we are defining a 
consistency volume around each burst trigger in channel H. If one of the mapped 
triggers from channel X falls in to this volume, we veto the H trigger. This is 
schematically illustrated in Figure 0] 



4-3. Software injections 

In order to estimate the variance of the distribution of HACR errors in estimating the 
parameters of the waveform, sine-Gaussian waveforms are injected into white noise 
and the HACR ETG is run over the data. The parameters estimated by the ETG are 
compared to the injected parameters in order to calculate the errors in the parameter- 
estimation. Distributions of the errors (after subtracting the mean and normalizing 
by the standard deviation) are plotted in Figure [5] along with the probability density 
of the Normal distribution with mean and variance 1. The standard deviation of 
the errors in the estimation of different parameters are plotted as a function of the 
SNR of the triggers in Figure It can be seen that, to a very good approximation, 
the standard deviation is an exponentially decreasing function of the SNR. The eight 
data points in the plots in Figure correspond to eight different values of the RSS 
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Rejection probability Rejection probability 

Figure 8. [Left]: Estimated veto efficiency plotted against the rejection 
probability. [Right]: Estimated false- veto probability plotted against the rejection 
probability. 

amplitude used for the injections. We fit the data with a simple polynomial fit and 
take it as the functional form of the standard deviation. 

We simulate two populations of bursts in channel X and if, as described in 
Section 13.11 The two data streams are processed by HACR and two sets of triggers 
{Ex} and {Eh} are generated. The X triggers are mapped to channel H using the 
transfer function from X to H . Figure shows the characteristic amplitude E a of the 
three population of triggers (Ex, E x and Eh) plotted against the central frequency 
. The injected events in channel X span 9 different amplitudes. The E x triggers 
and Eh triggers can be seen to fall nicely in to the shape of the transfer function. 

The veto analysis is repeated with different thresholds. The estimated veto 
efficiency is plotted against the rejection probability corresponding to the chosen 
thresholds in Figurc|Hl Although the estimated rejection efficiency roughly agrees with 
the predicted rejection probability, the effect of relying on a number of assumptions 
can be immediately seen. The main source of error in the analysis comes from the 
inaccurate parametrization of the errors in ETG in the parameter-estimation of the 
burst waveforms. This can be different for different ETGs. 

As described in Section im we also estimate the false- veto probability by injecting 
sine-Gaussian waveforms with random parameters in to the two data streams and 
preforming the veto analysis. The estimated false-veto probability is plotted in 
Figure |H| (right) as a function of the rejection probability. This exercise suggests 
that a veto efficiency of ~ 70% can be achieved at the cost of a false- veto probability 
~ 1%. 

5. Summary and outlook 

In this paper we have formulated and demonstrated a strategy to veto spurious noise 
transients that appear in the output of interferometric GW detectors. This method is 
novel in the sense that it makes use of the measured coupling between an instrumental 
sub-system and the detector output in order to veto a spurious trigger. The basic idea 
is that the noise recorded in an instrumental channel X can be mapped in to the 
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detector output (channel H) using the transfer function from X to H . This enables 
us to project the noise at the detector output in to two orthogonal directions - along, 
and orthogonal to, the mapped noise from channel X. If a non-stationarity in H is 
originated from channel X, it leaves the statistics of the component of H orthogonal 
to the 'mapped' Aunchanged. This can be verified by a statistical hypothesis testing. 
We have also proposed a less-rigorous, but computationally inexpensive alternative to 
the above method. In this method, the parameters of the triggers in channel X are 
compared to the parameters of the triggers in channel H to see whether a trigger in 
channel X can explain a trigger in channel H. 

Work is ongoing to incorporate the veto in the online data characterization 
pipeline of the GEO 600 detector [55]. So far, we have not explored the effect of the 
errors in measuring the transfer function. This certainly needs to be considered when 
the veto is applied to the actual data characterization/analysis pipeline. We leave 
that as future work. Also, the assumption that the transfer function is time-invariant 
is also not strictly true. Transfer functions in actual detectors can vary in time. This 
issue can be addressed by making repeated measurements of the transfer function, and 
tracking its evolution by continuously injecting and measuring spectral lines at certain 
frequencies. This is described in ^H] and is already being practiced in GEO 600. It 
may also be noted that the noise in the present-generation interferometers are not 
stationary Gaussian, and exhibits tails in the distribution. An actual implementation 
of this method should take this also in to account. 

The 'trigger mapping' veto needs to make certain assumptions about the power 
spectrum of the glitch in channel X. The assumption that we made in the paper, that 
the power spectrum can be approximated by a Gaussian function, should be verified 
against real-life glitches. It might also be possible to make assumptions which are 
closer to the reality, using better parametrization of the underlying waveforms. 

The results of the two veto methods presented in this paper may not be compared 
directly. The very choice of sine-Gaussian waveforms for the software injections 
especially favour the trigger mapping method, as they satisfy the assumption about 
the power spectrum of the burst waveforms that we used in the analysis. Yet, the 
results show that the 'noise projection' method performs better. It may be noted that 
the trigger mapping method also depends on the accuracy of the parameter-estimation 
of the ETG while the 'noise projection' veto relies only on the accuracy with which 
the transfer function is known. 
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